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C^ ■ I outline the effective field theory (EFT) calculation of nucleon-nucleon scatter- 

r^ ' ing which was recently carried out to next-to-next-to-leading order (NNLO) by 

Noam Shoresh and myself. In this calculation only potential pion contributions 
are included. These are the leading contributions from pions. A toy model was 
also considered in which potential pions are the only pion effects. The effective 
range expansion (ERE), which is valid for low energies, is used to derive matching 
conditions on some of the EFT couplings. Renormalization group flow analysis is 
used to fix the rest of the couplings leaving no free parameters at NNLO. 
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1 Introduction 

Effective Field Theory (EFT) is useful for systems with a clear separation of 
scales. This is the case in nucleon-nucleon scattering, where the light pions have 
q\ . a mass of m^ ~ 140MeV, and the next particle to consider, the p meson, only 

comes in at m p ~ 770MeV. In describing the two-nucleon system, however, 
special care is needed because of the appearance of a new low energy scale, 
namely the singlet channel scattering length. The Kaplan-Savage- Wise (KSW) 
power counting schemed was formulated to address this issue. 

In the EFT where the pions are included explicitly, the contribution of 
loops containing pions can be separated into (a) potential pions and (b) ra- 
diation and soft pions. Potential pions, corresponding to instantaneous pion 
exchanges, arise when the nucleons in the loop are put on-shell. Since both 
jj ' nucleons can be put simultaneously on-shell these contributions are enhanced. 

The radiation and soft pions enter when the pions are put on shell and describe 
retardation pion effects. 

In this talk I describe the contribution of the potential pions aloneH. It is 
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argued below, Subsection 8.2, that this is sufficient for a next-to-next-leading 
order (NNLO) calculation of nucleon-nucleon 1 S'o scattering amplitude at ex- 
ternal energies close to the pion production threshold. At energies much below 
this threshold, there could be contributions from radiation and soft pions. 
Nevertheless, the calculation with only potential pions serves as a probe for 
a better understanding of the KSW power counting, the appropriate fitting 
procedure and the role of renormalization group flow analysis. It is worth- 
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while to consider a model in which only potential pions exist, and one does 
not have to worry about the effects of the unaccounted for radiation and soft 
pions. Such is the two- Yukawa model described in Section g. The effective 
Lagrangian with potential pions for center of mass momentum p < m p /2 is 
described in Section |3|. KSW power counting is briefly reviewed in Section |f| 
The next-to-leading order (NLO) and NNLO scatteri ng a mplitude is calcu- 
lated in Section || and || respectively. In Subsection 7T , some of the EFT 



couplings are determined from low energy matching conditions for the ampli- 



tudes. Renormalization group (RG) flow analysis in Subsection 7.2 fixes the 



rest of the couplings. Results for the toy model arc presented in Subsection 3.1 



whereas in Subsection 8.2 I have included the corresponding results for real 1 5'o 
scattering. 

There are also relativistic corrections to the amplitude at NNLO due to 
the finite mass of the nucleons. These effects are very small compared to the 
potential pion contributions and can be neglected. 



2 The Toy Model 

We consider a toy model in which non-relativistic nucleons interact via two 
instantaneous Yukawa potentials: 



— m^r —m p r 

(r is the distance between the nucleons). This model is complex enough in 
the sense that the interactions involve two well separated scales. The Yukawa 
potentials describe the instantaneous exchange o£the n and a scalar p meson 
(without the complications of retardation effects cl). 

The realistic 1 Sq one pion exchange amplitude (here given in momentum 
space) is 
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where the first term is a contact interaction, and the second is the Yukawa 
part. The contact interaction is not included in the toy model, g^ is chosen to 
reproduce the Yukawa part. Comparing Eq. (0) and Eq. (pi), we set 
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The nucleon mass is taken to be M — 940 MeV, as well as m^ = 140 MeV, 
= 770 MeV, /„ = 132 MeV and g A = 1-25. g p is tuned to give a large 



scattering length (g p = 13.5 



-24.7 fm) by numerically solving the 



Schroedinger equation with this potential. 

The potential pion contribution to the field theory amplitude in this toy 
model and the real ^o nucleon- nucleon scattering are the same. One might 
wonder about the contact piece in Eq. (||) that was neglected. As will be 
explained in Section H, in the field theory amplitude the contact piece of the 
one-pion-exchange (OPE) can be absorbed in the definition of one of the 4- 
nuclcon couplings. 



3 The Effective Lagrangian 

The effective Lagrangian for the toy model appropriate for center of mass 
energies much smaller than m p is presented here. The effects of the short 
range force, p exchanges, are reproduced in the EFT by contact interactions. 
It is convenient Q to separate the effective Lagrangian in terms of the number 
of nucleon fields present: 

where C n describes the propagation or scattering of n nucleons. L\ is the usual 
kinetic energy term 

C 1 =N\id +V 2 /2M)N : 

where the isodoublet field N represents the nucleons: 
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Next, we write C 2 , the part of £2,-jthat has non-vanishing matrix elements 
between 1 Sq states. The matrices PkU are used to project onto the 1 Sq state, 
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where the a matrices act on the nucleoid spin space and the r matrices act on 
the nucleon isospin space. The result isn't 
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(summation over the repeated isospin index i is implied.) 

We use dimensional regularization where the Lagrangian is given in D 
space-time dimensions and /x/2 is an arbitrary mass scale introduced so that 
the couplings C2n, -D^n have the same units in any dimension D. Cin are 
the coefficients of contact interactions containing 2n derivatives and -D^n 
are coefficients of contact interactions involving 2m derivatives and 2(n — m) 
powers of the pion mass. We use the convention D 2n = &%„, Ci n = D^n ■ 
The ellipses indicate operators with higher powers of derivatives or pion mass 
insertions. 

The interaction with the pions is described by a non-local term in the 
action: 

d D x d D y (N T (x)P i N(x)) T g OT * _ 5(x° - y°) (N T (y)PN(y)) 

d D p d D P < (N T (p)P t N( P ))' A ^ (N^p^PNip')) . (6) 

Once the effective Lagrangian has been described, one needs a power count- 
ing scheme, i.e. a set of rules that determine the relative sizes of diagrams 
contributing to the scattering process. 

4 The power counting 

KSW power counting is motivated by the large scattering length a. It can 
be better understood by looking at the singlet channel amplitude at very low 
energy where even the pions can be treated as heavy particles and hence in- 
tegrated out. The effective field theory amplitude is related to the phase shift 
by 

A= ^L l , (7) 

M p cot 8 — ip 1 



where p = y/ME is the center of mass momentum and S is the singlet channel 
phase shift. For low energies one can write U the effective range expansion 

(ERE): 
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Here the scattering length a can be arbitrarily large, while the other co- 
efficients ro, ri, ... are assumed to be of natural sizes, r n ~ 1/A 2n+1 , where 
A ~ m n is the cut-off for this low energy theory, p is taken to be small 
compared to the high energy cut-off A but not compared to 1/a <C A. The co- 
efficients in the two different ERE expansions are related perturbatively. The 
first two terms are 
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The difference between 7 and 1/a enters at NLO, whereas the difference be- 
tween s arL d Tq is only important at orders higher than NNLO. 

Substituting Eq. (g) in Eq. (Q), it can be seen that the amplitude has a 
pole at p = ij and this sets the radius of convergence for the Taylor expansion 
about p = 0. Formally identifying 7, p ~ Q, where Q/K is the expansion 
parameter, the amplitude Eq. (M) can be expanded about the pole in powers 
of p/K and 7/A, 
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In KSW power counting the field theory amplitude is organized as an ex- 
pansion in Q, Eq. (UQ). It has been shownQ that the leading term in expansion 
Eq. (110) is obtained by summing up the bubble chains in Fig. Fy with insertions 
of Co, giving the leading order (LO) amplitude: 
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Figure 1: Leading order contribution A— 1 from non-perturbative insertions of Co- 

where L is the loop integral 
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A subtraction scheme that reproduces the KSW power counting is the 
power divergence subtraction (PDS) scheme that was introduced mo. In PDS 
one subtracts the poles of L in both D = 4 and D = 3 dimension (which gives 
X su b = \x in Eq. (|l2|)). Setting /i ~ p, all the diagrams contributing to the 
bubble chains are rendered equal in size (~ 1/p)- This justifies summing up 
the bubble graphs at LO. 

In the theory with pions, pion effects are assumed to be perturbative. From 
the PDS scheme follows a power counting scheme - KSW - in which: 

• The expansion parameter is Q/A, where p, /i, m^ ~ Q. 

• >s\. 2-jn=— i n ' n ^ ^ ■ 

• The renormalized coupling Co = Co(/u) scales as l//i ~ 1/Q, and more 
generally the renormalized couplings C2n, D%™' scale as 1/Q n+1 . 

Some immediate consequences are: 

• The loop integral L scales like p or //, i.e. O(Q). Therefore, adding a C 
and a loop L to a diagram does not change its order since the powers of 
Q cancel. It is now required to dress each diagram in the theory by Cq to 
all orders. (This is diagrammatically represented by the "blob", Fig. g). 

• When the nucleons are put on shell, nucleon propagators scale as 1/Q 2 



and loop integrals J dq A scale as Q 
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Figure 2: Definition of the "blob" . 

This theory should correctly describe the physics at momentum below 
A ~ rrip/2, which is the scale associated with the cut in the amplitude due to p 
exchange. A calculation to order Q n is expected to be accurate up to an error 

or~(Q/A)" +1 . 

It is a key feature of the KSW power counting that the LO amplitude 
has the right pole, or the correct scattering length. However, at higher orders 
there are contributions that move the position of the pole in the amplitude 
and counterterms must be introduced to cancel, these effects. This is done by 
expanding the EFT couplings in powers of QellM: 

Co = Co,_i + Co,o + Co.i + • • • 

C% = C<2 } —2 + C/2,— 1 + • • • 

Ci = 04,-3 + C*4,-2 + ' ■ ■ 

; (13) 

where the second subscript denotes the scaling with powers of Q. 

In the last three sections, all the components of the computational frame- 
work of the EFT for nucleons have been identified: the effective Lagrangian 
(Feynman rules), the regularization and subtraction methods, and the power 
counting scheme that organizes the diagrams of the theory in a perturbative 
series. In the following sections we apply this formalism to the calculation of 
the So scattering amplitude. 



5 The Amplitude at NLO 

In the theory with truly dynamical pions, the contact piece of the pion exchange 
can be re-absorbed in the definition of Co,o- Hence, there is no difference in the 
potential pion amplitude for this toy model and the real N-N scattering am- 
plitude. The potential pion contribution at NLO has already been calculated 
beforetx I include the calculation here for completeness and clarity. 

At NLO, we need to consider local operators that make up the following 
four-nucleon vertices, classified here according to their Q counting: 

OiQ- 1 ) : -Co.-x 

O(Q0) : _(s!+sHc 2) _ 2 ;-C ,o;-rr4D 2 ,_ 2 

with qi , q 2 the single nucleon incoming and outgoing momenta, in the center 
of mass coordinates. Only the linear combination Co,o + tiJ-C>2,-2 appears 
in the amplitude, and it is denoted by -Bo.o- From this point on we adopt 
the convention that the second subscript is omitted from leading couplings, so 
that, for instance, Co stands for Co,-i, C 2 is really 6*2,-2, etc. 

The diagrams contributing to the NLO amplitude are shown in Fig. pi 
The amplitude is 

where A includes all the diagrams with no pion exchange, and -4 con- 
tains all the diagrams with a single pion exchange. It is convenient to express 
Fig. H(b) and Fig. BKc) in terms of the loop integrals P\ and P 2 , so that 
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Figure 3: Next-to- leading order contribution Aq. 
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Eq. (ITq) , where only the poles have been subtracted, differs from an earlier 
calculation of P2U that involved finite subtractions. 

6 The Amplitude at NNLO 

At this order, there are insertions of three kinds of local operators, listed here 
according to their Q counting: 
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Figure 4: Next-to-next-to-leading order contribution A^ from diagrams with no pion 
exchange. 



©(Q" 1 ) : -Co 

O(Q ) : -Bq.o;- 15 ^*?: 



©(Q 1 ) : -B 0)1 ;-M+^lB 2i _ i; _p*C4 

There are 3 new couplings that enter at this order: i?o.i = Co,i + m%D2,-i 
m*D 4 , B 2 .-i = C 2 _i + m^-Of ) and C 4 . The NNLO amplitude is 

A x = A { t ] + A^ + A { r ] ■ 
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6.1 A\ 

The diagrams in Fig. are generated by two insertions of 0(Q°) operators or 
a single insertion of 0(Q l ) operator, and give: 
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ff.jg -A^ 

All the diagrams of Fig. are generated by insertions of local operators of 
order 0(Q ) and a single pion exchange. The one pion exchange contribution 
at this order is 
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Figure 5: Next-to-next-to-leading order contribution .4^ from one pion exchange diagrams 
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Figure 6: Next-to-next-to-leading order contribution A.\ from two pion exchange dia- 
grams. 



and Pi t 2 are defined in Section @. In carrying out the subtractions in Qi and 
Q2, some care is needed. The poles and the corresponding counterterms should 
be calculated treating Co perturbatively. 



6.3 4 7 " 

Fig. rcj contains the diagrams of this set. The graphs Fig. pKe), (d) and (f) can 
be written in terms of analytic functions Ii(p), Iz(j>) and I^ijp) respectively 
which are presented in the Appendix. The rest of the diagrams can be written 
in terms of expressions that have already been defined. We find 
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In the last two sections, we have calculated an expression for the LO+NLO 
+NNLO amplitude in terms of six independent free parameters Co, -Bo.O: C2, 
Bq.I: -62,-1 and C4. The table below shows the new free parameters introduced 
at each order in the expansion. 



LO: 


Co 






NLO: 


-Bo.o 


c 2 




NNLO: 


B0.1 


-62.-1 


c 4 



When the number of parameters is doubled (from NLO to NNLO) it is hard 
to consider an improvement of the fit as evidence for the convergence of the 
EFT expansion. Therefore, it is desirable to find theoretically viable ways to 
reduce the number of independent parameters. This is the subject of the next 
section. 

7 Determining the Effective Field Theory Couplings 

7.1 Matching conditions 

In this subsection I discuss a method for deriving conditions that reduce the 
number of free parameters. The guiding principle is simple: the leading cou- 
plings are defined in such a way that certain features of the low energy phase 
shift (like the effective range expansion parameters) are reproduced exactly. 
Conditions are then imposed on the sub-leading couplings so that these fea- 
tures are not jnodified at higher orders. Similar matching conditions were 
introduced inQ. 

From Eq. (J7|) , the Q expansion of p cot S in terms of the amplitude A — 
A-i +A +Ai-\ , is 
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At low energies, the same quantity is given by expansion Eq. (g 
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where pcotS has been expanded about p 2 = — 7 2 (or p = ±17). 

The LO amplitude contributes only a constant to p cot S. Cq is determined 
by requiring 

47T 
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This puts the pole in the amplitude at p = vy. At NLO and NNLO there are 
more constant contributions to pcotd. The requirements that the position of 
the pole does not shift at NLO, nor at NNLO, reduce to: 

NLO : A|p=-* 7 = 
NNLO : -4i|p=_iT = 0. (25) 

These conditions determine the couplings _B ,o an d ^o,i- Similarly, one requires 
that so in Eq. ( f23|) is determined at NLO by C2: 

8 ( At: A \ 1 

V {-ma^aTi) [p2 =-~< 2 = 2 s0 ' (26) 

and that higher order (NNLO) effects do not change s^: 

NNLO: |-A|p=-n=0. (27) 



This condition determines B2.-1. (Eqns. (|25[) and ( |27| ) can be more directly 
obtained by considering the expansion of the amplitude Eq. (|10|).) Note that 
at NNLO, fixing 7 is not exactly the same as fitting the scattering length a 
whereas reproducing sq is essentially identical to fixing the effective range rg . 
The values of 7 and sq are used to determine Co and Ci. One can determine 
C4 from si, the next frea-Barameter in the expansion Eq. (123) . However, as in 
the theory without pionsEJo, it can be shown using RG analysis that at NNLO, 
C4 is completely determined. This is done below. 

7.2 Renormalization group flow 

The RG equations are derived by requiring that the amplitude be exactly \x 
independent at each order in Q for any value of p. The solutions to the RG 
equations for Co and C2 are: 

4tt / 1 x ! 



C 3W = fo («) 2 . (28, 

The integration constants £c and £c" 2 are determined by boundary con- 
ditions. The other coupling constants, denoted here by X, have the general 
solution 

*(A0=fr(^§^) +/x(Co(A*), J Bo,o(^),C 2 ( M ); Ai ) (29) 
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where fx is known, and the only freedom in X is in the constant ^x, deter- 
mined by boundary conditions. Specifically, 

2B Q . C 2 fMCaV 

Examination of the solutions in Eq. ( |30| ) shows that the \i scaling of some 
couplings do not agree with their expected Q counting. To determine the 
Q counting one also needs to know the Q counting of the RG constants £x- 
That 7 can make up Q counting for some of the £xs was already shown in 
the theory without pions, where, for example, C2.-1 ~ 7 2 /(/ i — 7) 3 B'Q. It 
is possible, however, to determine the Q scaling of the leading couplings. We 
demonstrate this point by discussing directly the coupling C4, which has power 
counting of 1/Q 3 . The second term in the solution to the RG equation for C4 
(Eq. (poj)) is ~ £c 4 /m 2 - I n order for it to contribute at NNLO, £c 4 should be 
0(1/ Q). Since C4 has a sensible large scattering length (7 — ► 0) limit, the 
possibility £c t ~ 1/7 is ruled out. It is also assumed that C4 is associated 
with short distance physics and thus it is fundamentally independent of pion 
physics. These considerations allow setting £c 4 = at this order, and have C4 
be completely determined by Cq and C2 . There are in fact no new parameters 
at NNLO! 

In the next section we check how well the EFT reproduces the phase shift 
of the toy model and the real data when the matching conditions derived in 
the last two subsections are imposed. 

8 Results 

The 1 Sq scattering cross section a is related to the phase shift and amplitude 
by 

4tt . 2 
(7 = — sin 

M 2 
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which can be expanded in powers of Q ( writing a = er_2 + o~—\ + o"o + ' ' ') : 
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M 2 
ff ° = ^ (-A-i^i + X;A + | A| 2 ) • (32) 

There are several advantages to using a as opposed to 5 in comparing EFT 
result to data. The experimentally measured cross section a can be related to 
the sizes of the absolute value of the amplitude Eq. (p^). 5 is very sensitive to 
the value of n and other higher moment coefficients and it is reflected in the 
value of Ci that one chooses. This problem is avoided if one plots a or sin 2 5 
instead, where the effects of n or £c 4 are truly sub-leading. 

In the following subsections we compare the EFT result with the exact 
phase shift from the toy model, and also do the corresponding comparisons for 
the real 1 S phase shift. We consider how well the data is reproduced by the 
EFT, whether the expansion is converging, and what is the breakdown scale. 



8. 1 Results for the toy model 

For the toy model, 7 = —7.743 MeV and so — 1.584 fm and these values are 
used to determine the couplings. The results are presented in Fig [?]. Fig [7] (a) 
shows a clear improvement of fit as we go to higher order in the expansion. The 
NNLO result is a prediction-as we do not introduce any new free parameter. 
Fig. [7] (b). shows the errorsEa at each order, A = | sin 2 (5) eft — sin 2 ((5) ea;act |, 
on a logdog scale. The breakdown scale, where the errors become comparable, 
can be read off of the error plot to be ~ 520 MeV which is slightly higher than 
the expected m p /2. 

Table [H shows the couplings at n = m^. The entries are given in units 
of fm . To make a direct comparison of the couplings easier, we included 
the momentum factors from the operators C2P 2 , B2.-1P 2 and C4P 4 , evaluated 
at p = m-x- From the table we see a hierarchy in the sizes of the couplings 
corresponding to the orders in the expansion. This is a verification of the power 
counting scheme. A rough estimate of the expansion parameter is Q/A ~ 
rriTr/imp/2) ~ 1/3, which is also suggested by the table. 



1G 



p (MeV) 



Figure 7: Toy model: (a) sin 2 (<5) vs. p(MeV); (b) log-log plot of A e sin 2 (8) E FT - 
sin 2 (6) exa ct\ vs. p(MeV). The exact phase shift is described by the solid curve. The long- 
dashed curves denote the LO phase shift and error, the dashed curves are the NLO results, 
and the dot-dashed curves describe the NNLO results. 



Table 1: Renormalized couplings for the toy model at fj, = m^ in units of fm . 



LO 


NLO 


NNLO 


Co 


C 2 ml 


Bq,o 


Bq,i 


£2,-1 m'i 


Ca m% 


-3.526 


0.06912 


1.412 


0.5620 


0.2886 


-0.001355 



8.2 Results for the real 1 Sq scattering 

Let us recall that in the KSW power counting the LO amplitude, which is 
inversely proportional to the momentum p, is calculated in terms of a non- 
perturbative dimension 6 four-nucleon operator. The other higher derivative 
operators and pion exchanges are treated perturbatively. The expansion pa- 
rameter is Ql h. with p, TO?,- ~ Q and A some high energy cut-off, and at 
p ~ Q the LO amplitude scales as 1/Q- Pion production, however, at exter- 
nal momentum p — ^/m^M ~ 360 MeV, introduces a new scale, and it is 
convenient to use^ different expansion parameter, Q r ~ y/m^M , instead of 
Q ~ Ql/M ^C Q r Q- At p ~ Q r , the LO amplitude scales as 1/Q r . The relative 
importance of the potential, radiation and soft pion contributions depends on 
the physics associated with different scales determined by the momentum p. 

At p ~ Q r , potential and soft pions contribute at Qjl and Q^. respectively, 
while radiation pions contribute at Qf H Thus, for an NNLO (0(Q r )) cal- 
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Table 2: Renormalized couplings for the real 1 Sq channel at fi = m^ in units of fm 



LO 


NLO 


NNLO 


Co 


C 2 mi 


-Bo,o 


-Bo,i 


-82,-1 mi 


C 4 m% 


-3.5226 


1.39685 


1.41738 


0.56156 


0.541789 


-0.553904 



dilation, only potential pions contribute at this scale. At p ~ Q, potential 
pions contribute at order Q°. The power counting for the soft and radiation 
pion is less clear. Examination of the radiation pion diagrams that have been 
calculated in El, shows that some diagrams, which scale as Q^. at high mo- 
mentum, become order Q 1 ' 2 at low external momentum. Such contributions 
seem to violate the presumed expansion of the amplitude in integer powers 
of Q, but in fact there is a cancellation between the different diagrams, and 
the leading contribution from radiation pion graphs starts at 0(Q). Further- 
more, the remaining terms are highly suppressed because of Wigner's SU(4) 
spin-isospin symmetry a and can be neglected in the NNLO calculation (see 
Tom Mehen's discussion in these proceedings). We do not see (at the moment) 
a general reason for such cancellations (of order Q 1 ^ 2 ), and it is in principal 
possible that radiation or soft pion graphs that are of high order in Q r at 
high momentum, become important at low momentum, and should in fact be 
included in an NNLO calculation of the amplitude. Here, however, this option 
will not be investigated, and only the potential pion contributions will be taken 
into account, admittedly with the possibility of having left out some graphs. 
The EFT amplitude for the toy model with potential pions is identical to the 

potential pion contribution to the real 1 Sq amplitude, after a redefinition of 

2 
Co - * Co — nrz- With these considerations, the EFT amplitude with only 

potential pions is applied to the real 1 S'o scattering. We use 7 = —7.889 MeV 
(a = —23.68 fm) and sq = 2.73 fm to fix the couplings. 

Here again, we have a prediction at NNLO, Fig. g. By going to higher 
order we get better result, Fig. |](a). From Fig. ||(b), the breakdown scale is 
estimated to be ~ 400MeV. Just like in the toy model, a rough estimate of the 
breakdown scale is m p /2 due to the exchange of vector meson p. 

Table shows the couplings at /i = m^. The entries are given in units of 
fm . We included the momentum factors from the operators C2P 2 , i?2,-ip 2 



and C4P , evaluated at p •■ 



The sizes of the couplings at various order in 



the expansion suggest convergence of the expansion with Q/A ~ 1/3. 
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Figure 8: Nijmegen: (a) sin 2 (<5) vs. p(McV); (b) log-log plot of A 5 \sin 2 (8)EFT — 
sin 2 (<5)jvijme9en I vs - p(MeV). The Nijmegen phase shift is described by the solid curve. 
The long-dashed curves denote the LO phase shift and error, the dashed curves are the NLO 
results, and the dot-dashed curves describe the NNLO results. 



9 Summary and Conclusion 

A NNLO calculation for nucleon-nucleon scattering in the 1 5'o channel using 
potential pions was presented. Radiation and soft pions do not contribute at 
this order for center of mass momentum p ~ Q r . However, there could be 
soft and radiation pion contributions at momentum p ~ Q which have not 
been calculated yet. If for low momentum p ~ Q, soft and radiation pion 
contributions are found to be small at NNLO, then this calculation gives the 
NNLO nucleon-nucleon scattering amplitude in the Sq channel for momentum 
P < m p /2. 

The calculation was applied to a toy model and the real data. The toy 
model provided a setup where only potential pions contribute by construction. 
This allowed us to address some general questions regarding the convergence 
of KSW power counting, number of free parameters and the role of renormal- 
ization group flow, without worrying about unaccounted for radiation and soft 
pion contributions. 

At NNLO there are six independent parameters. In the effective field 
theory, these parameters are associated with contact operators which encodes 
information about the high energy physics. However, in the singlet channel 
the high energy physics conspire to produce a low energy scale 7^8 MeV. 
This small energy scale has to be taken into account in using dimensional 

19 



analysis to determine the sizes of the parameter. RG analysis, which only 
determines the fj, scaling of parameters, is not always enough to determine 
their Q scaling. Some of the parameters such as £c 4 can be determined from 
RG analysis. We introduce a method to determine the rest of the parameters. 
This is done by matching the EFT amplitude at low momentum to the effective 
range expansion and all the parameters are given in terms of two experimental 
inputs, 7 and s - We find that there are no free parameters at NNLO and we 
make a prediction for the high momentum behavior of the scattering amplitude. 

For the toy model we find that including higher order terms in the ampli- 
tude improves the fit to "data" . The hierarchy in the sizes of couplings show 
converges with an expansion parameter of Q/A ~ 1/3. The breakdown scale is 
estimated to be ~ 520 MeV by comparing the sizes of errors at different order. 
This is slightly higher than the expected scale m p /2. 

The fit to real data shows significant improvement at NNLO. We see an 
hierarchy in the sizes of the couplings which suggests an expansion parameter 
similar to the toy model Q/A ~ 1/3. The breakdown scale is estimated to be 
~ 400 MeV. 
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Appendix: The integrals I\, I<x and I3 

Fig. @ (e) is given by: 



<f- 



2 



h = 



d 3 q 1 1 1 



(2-7r) 3 q 2 — p 2 — ie (q — p) 2 + m 2 (q — p') 2 + m 2 

(iL(p) + M{p)) (33) 



327rm 3 /} 3 
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L(,)= 1 L^i^m A1 + Vt)(1 + V(l " /,) 



t(i - 1) \ 1 + v 



iln 2 (l + 4p 2 ) 



1/2 

M(p) = 2 / tft \ (ta,n- 1 (2p^l + 4 P H(l-t)) - tan" 1 ^)) 
o 

= - 7T 2 - z In 2 (2) + - In (4) In (1 - 2 i p) + 2 ?r In (2 - 2 i p) 
-i In (2 (1 + i p) 2 (1 - 2 i p)) In (1 + 2 i p) 

- tan' 1 (2 p) In f - |^g J - 2 tan" 1 (p) In (l + 4 p 2 ) 

-2iii 2 ( -+ip) -2iLi 2 (-l-2ip)-2iLi 2 (2-2ip), 

where the average is over the outgoing momentum p'. p = p/m^ and the 
dilogarithm, Li 2 {z), is defined as 



Lh{z) = -ldt Hl -*\ 



Fig. H (d) is given by: 



2 \ 2 
9 \ _4,, f 2; 



d 3 q d 3 l 1 



(34) 



h = 



(27r) 3 (27r) 3 q 2 -p 2 - ie (q - p) 2 + m 2 I 2 - p 2 - ie (q - 1) 



2 



dxln In 



128n 3 m 2 p J_ QO \i — x + p J \ (x — p) 2 + \) x 2 — p 2 — ie 

__J-^(-^ + 4z^ln(2-2zp) + ln(2)ln(l + 2*p) 

+ In (1 - 2 i p) [4 In (1 - i p) - In ((1 - 2 i p) (l + 4 p 2 ) /8)] 
-Lz 2 [ 2~ ip ) +Li2 ( 9 +ip ) + 4L *2(2-2ip) 
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From Fig. p](f), we define 



2 \ 2 
9~ \ 4 



'■'■■ = ^ \Sp) <m3!s 

d 3 k 1 

J2nf k 2 -p 2 - ie 
d 3 q 1 



h = 
/(k,p) = 



/(k,p) 2 
1 



(2ir) 3 q 2 — p 2 — ie (q — k) 2 + m 2 



and it gives 



1/2 



64ir 3 p 



dt 



1 



t(l-t) 



dse -(m- 2l p) S 



hi 



1 — ip 



647r 3 pm 7r \_ ^1 — «2/9 
+Li 2 (l - ip) - Li 2 (l - i2p) 



ln(ip) + Li 2 



-i2pst 1 



-IP 

l-i2p 



(35) 



(36) 



The explicit expressions for I\, I 2 and I3 presented above are valid only 
for p > and small positive imaginary p (\p\ <C 1). For other values of p, 
the integrals have to be done taking into account the poles and cuts in the 
integrand appropriately. 
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